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Abstract 

We analyze quantum quenches in integrable models and in particular the determination of the 
initial state in the basis of eigenstates of the post-quench hamiltonian. This leads us to consider 
the set of transformations of creation and annihilation operators that respect the Zamolodchikov- 
Faddeev algebra satisfied by integrable models. We establish that the Bogoliubov transformations 
hold only in the case of quantum quenches in free theories. In the most general case of interacting 
theories, we identify two classes of transformations. The first class induces a change in the S- 
matrix of the theory but not of its ground state, whereas the second class results in a "dressing" 
of the operators. As examples of our approach we consider the transformations associated with a 
change of the interaction in the Sinh-Gordon and the Lieb-Liniger model. 
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1 Introduction 



A quantum quench is an instantaneous change in the parameters that determine the dynamics of an 
isolated quantum system e.g. the masses or couphng constants of its hamiltonian. This topic has 
recently attracted a lot of attention as shown by the increasing number of papers addressing this 
issue (for a recent review, see [T] and references therein) . From an experimental point of view this 
is a feasible way to bring the system out-of-equilibrium and study its evolution under the quantum 
mechanical natural laws, in isolation from the environment. In particular, the scientific interest 
in quantum quenches started growing after the experimental realization of global sudden changes 
of the interaction in cold atom systems, a novel technology where quantum statistical physics can 
be experimentally demonstrated and probed [2H5]. From a theoretical point of view the problem 
consists in preparing the system in a particular trial state, which is typically the ground state of 
some hamiltonian, and study its evolution under a different hamiltonian [Sl417| . Apart from being 
one of the simplest and well-posed ways to study out-of-equilibrium quantum physics, quantum 
quenches also give rise to a fundamental long-standing open question of central importance in 
statistical physics, the question of thermalization: how do extended quantum physical systems 
tend to thermal equilibrium starting from an arbitrary initial state? 

Of particular interest is the case of (l+l)-dimensions where a discrimination between integrable 
and non-integrable systems is possible. Integrable models are models that exhibit factorization 
of the scattering matrix and can be solved exactly (see, for instance, [18] and references therein). 
Their classical counterparts possess as many integrals of motion as their degrees of freedom and 
this fact prevents thermalization of an arbitrary initial state, as not all of the micro-states of equal 
energy respect the conservation of all other integrals of motion. This property is also expected 
to hold at the quantum level. In a seminal experiment [5] it was observed that a trapped (l+l)d 
Bose gas, initially prepared in a non-equilibrium state, does not thermalize but tends instead 
to a nonthermal momentum distribution. The absence of thermalization suggests as a possible 
reason the integrability of the system which approximates a homogeneous (l+l)d Bose gas with 
point-like coUisional interactions, a typical integrable model, even though the confining potential 
used in the experiment breaks the homogeneity and therefore integrability of the system. This 
experiment triggered an intense discussion about the role of non-integrability in the thermalization 
process. It was soon conjectured that in an integrable case the system does exhibit stationary 
behavior for long times, described however not by the usual Gibbs ensemble but a Generalized 
Gibbs Ensemble (GGE) where new Lagrange multipliers are introduced into the density matrix, 
one for each integral of motion, for accounting their conservation (in the same way that the inverse 
temperature /3 is the Lagrange multiplier corresponding to the constraint of energy conservation) 



This conjecture has been shown to be correct in many different special cases, both by analyti- 
cal and numerical methods [TTldDHlS]- On the other hand, it has not become yet clear whether 
non-integrability alone is sufficient to ensure thermalization or not |27ll28j neither has exact ther- 
malization been firmly demonstrated as an outcome of unitary evolution. For instance, recent 
analysis suggests that the behavior may be more complicated and may depend on the initial state, 
finite size effects and locality |29H31| . For recent experimental developments backed by numerical 
simulations we refer to [32l434| . 

Concerning with the analytic approach to the problem, in the paper |24| it was shown that 
any quantum quench in an integrable quantum field theory where the initial state has the form 



leads to stationary behavior as described by the GGE ansatz. In the expression above, |0) is the 
vacuum state of the theory, while Z^{6) is the creation operator of the quasi-particle excitation, 
which satisfies the relativistic dispersion relations E ~ mcosh9, P — msinh^, where 9 labels the 
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rapidity of the particle of mass m. As it is clear from this expression, the function K{6) is the 
amplitude relative to the creation of a pair of particles with equal and opposite rapidity. 

The above form of quench state is also called "squeezed coherent" state. The reason to choose 
such an initial state comes from its relation with boundary integrable states (i.e. boundary states 
that respect the integrability of the bulk theory) and from the technical advantages it exhibits. 
It is however true that this requirement is satisfied in general for quantum quenches in a free 
theory, bosonic or fermionic, as well as for the important cases of Dirichlet and Neumann states in 
integrable field theory. These states are supposed to capture the universal behaviour of all quantum 
quenches in integrable models, if renormalization group theory expectations are also applicable 
out-of-equilibrium. They have been successfully used earlier pi1,ll2| . applying a Wick rotation 
from real to imaginary time which allows a mapping of the original quantum quench problem to 
an equilibrium boundary problem defined on a euclidean slab with boundary conditions both equal 
to the initial state right after the quench. Although this approach does not help in determining 
the expression of the initial state as a function of the quench parameters, it has led to correct 
predictions in certain important asymptotic limits. 

It still remains to find out from first principles, whether this assumption for the form of the 
initial state holds in general for any quantum quench in an integrable system or, if not. under 
what conditions this happens. Our method to attack this problem begins by understanding the 
fundamental reason why this condition holds generally for free systems and then investigating if 
this reason can be generalized to the integrable case. It turns out that in free systems the reason 
lies in the fact that the relation between the creation-annihilation operators before and after the 
quench is of linear Bogoliubov type, which itself is a consequence of their canonical commutation 
or anti-commutation relations. In integrable theories these commutation relations are replaced 
by the so-called Zamolodchikov-Faddeev (ZF) algebra which, assuming for simplicity that there is 
only one quasiparticle in the theory, can be written as 

z{ei)z^{92) = s{92 - 9i)z\e2)z{ei) + ^(^i - 62) (3) 

Intuitively this means that the exchange of two quasiparticles is done by the scattering matrix S{6). 
Then a natural question arises: what are the possible transformations of creation-annihilation 
operators that respect the above algebra? This is a question of more general interest in both 
the abstract mathematical description of integrable field theories and their potential physical 
applications in concrete models. In this article we show that, unlike in free theories, the ZF 
commutation relations do not admit Bogoliubov transformations and we construct several other 
classes of non-trivial infinitesimal transformations. 

We start our presentation by first discussing the structure of initial states in global quantum 
quenches. Then we outline a general strategy on how to determine the initial state from the re- 
lation between the ZF creation-annihilation operators before and after the quench and derive the 
conditions that must be satisfied by infinitesimal transformations of these operators in order to 
respect the ZF algebra. After showing that in interacting theories the linear Bogoliubov transfor- 
mations do not leave the ZF algebra invariant, we find two types of acceptable transformations, 
the first of which induces a shift in the ^-matrix but does not affect the ground state of the theory, 
while the second does not change the 5'-matrix but it does change the ground state and the cor- 
responding transformed ground state is, under certain conditions, of the squeezed coherent state 
form. We note that these are special classes of transformations and outline how more general ones 
can be constructed. Next we apply these ideas in two typical integrable models, the Sinh-Gordon 
model and the Lieb-Liniger model, deriving some examples of infinitesimal transformations of their 
ZF operators that demonstrate the presence of the first type constructed before. Finally we sum- 
marize our findings, giving directions for their application to concrete quantum quench problems. 
There are also two appendices: in Appendix A we discuss the squeezed states in free quantum 
field theories while in Appendix B we discuss the derivation of the two classes of generators of the 
ZF algebra transformations. 
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2 On the initial states in global quantum quenches 



A quench process consists of preparing the system in a state l^po) that is not an eigenstate of its 
hamiltonian H and let this state unitarily evolve according to H. At time t, the expectation values 
of local observables A(r) are given by 

(A(i,r)) = (V-ole^^^AWe-'^'IV-o) , (4) 

with similar expressions for higher point correlation functions. As evident from the expression 
above, important information on the subsequence dynamics of the system is encoded in the initial 
state \ipo). Relevant features of this state can be derived on the basis of general considerations 
for extended quantum systems having particle excitations and for global quenches. First of all, 
by relativistic (or even galilean) invariance, we can always assume that the quench state \iPq) 
carries no momentum. Let (p) be the creation operator of a particle excitatior0 of the system 
of momentum p and let us assume that a basis of the Hilbert space is given by the multi-particle 
excitations, eigenvectors of the hamiltonian H. Then the general form of the initial state Itpo) for 
a global quench is given by an infinite superposition of multi-particle states of zero momentum 

oo ^ n 
IV'o) = ^ dpi ... dpniCn{pi,-.-,Pn)S(^Pr)Z\pi)---Z'^{p„)\0) , (5) 

Tl = 0'' 1 = 1 

where |0) is the vacuum state of the system. This requirement is due to a thermodynamics 
argument related to the formulation of quench dynamics in d dimensions to the thermodynamics 
of a (d + 1) dimensional field theory in a slab geometry, where the initial state \iPq) plays the 
role of boundary conditions on both borders of the slab |12| . In this interpretation of the quench 
process, the quantity 

Zo(t) = (V-ole-^^lVo) = e-^"("' (6) 

plays the role of the partition function of the system with boundary conditions fixed by l-iAo)- 
For global quenches, the corresponding free energy J-oiT) must be an extensive quantity, J-o{t) ~ 
1^/o(t), where V is the volume of the system. On the other hand, this quantity can be computed 
by employing the expression ^ of the initial state \ipo) and a proper normalizatioijl of (5(0): then 
the only way to have an extensive behavior in the volume V of the system for J-o{t) is by jV'o) 
containing an inin7\nite number of multi-particle states. 

Notice that one way to automatically take into account the zero momentum condition of the 
initial state \tPq) is to assume that its infinite superposition is made of pairs of particles of equal 
and opposite momentum, i.e. "Cooper pairs" 

oo „ 

I'/'O) = ^ / rfpi . • • dpnJC2nipi, ■ ■ . ,Pn)Z\-pi)Z^ipi) ■ ■ ■ Z\-pn) Z\p„)\0) . (7) 
n=0 

It should be stressed, though, that this formula is a particular case of the more general form 
But even with this simplification, to specify the initial state j-iAo) one still needs an infinite number 
of amplitudes K,2n{pi, ■ ■ ■ ,Pn)- The great technical advantage of the squeezed coherent states, 
whose concise expression is given by 

l^o) ^ exp (^j dp K{p)Z\p)Z\-p)^ |0) , (8) 

becomes then evident. In this case, in fact, all the multi-particle amplitudes /C2n(pi, ■ • ■ iPn) can 
be expressed in terms of products of the single amplitude K{p) entering (|5]), therefore greatly 
simplifying the problem. 



^To simplify the following formulas we assume that the system has only one neutral elementary excitation, with 
normalization proportional to the <5 function, {Z{p\)Z'^ {p2)) oc S{pi — p2). 

^See 1351 for an explicit regularization of this term in the one-dimensional case. 
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Squeezed coherent states naturally appear in two contexts: (i) in the purely boundary inte- 
grable field theories considered by Ghoshal and Zamolodchikov [321, where the amplitude K{p) 
also satisfies additional conditions (boundary unitarity and crossing symmetry) and (ii) in quench 
processes in free theories, both bosonic and fermionic. In the latter case, it is worth noticing that 
the commutation or anti-commutation relations of the annihilation and creation operators Z{p) 
and Z^{p) of these theories can be cast in the form of ZF algebra ^ with = 1 for the boson 
and S = —1 for the fermion. The only parameter entering these theories is in this case the mass 
of the their excitation and, as shown in detail in Appendix A, its sudden change can be taken into 
account by a Bogoliubov transformation of the annihilation and creation operators. Since the Bo- 
goliubov transformations leave the commutation or anti-commutation relations invariant, in turn 
they can be seen as the transformations which leave invariant the ZF algebra of free theories. This 
observation leads us to investigate a more general class of transformations of the ZF operators in 
interacting integrable field theories which leave their algebra invariant. 

3 Quenches in Integrable Systems 

In this section we analyze the quantum quenches in systems which are integrable before and after 
the sudden change of one parameter Q, which can be for instance the mass of the particle or the 
coupling constant of the theory. One of the main issues of this problem is to write down the pre- 
quench state (usually the vacuum, annihilated by the pre-quench particle operators) in terms of the 
post-quench particle basis. This task involves in principle the computation of an infinite number of 
inner products, an operation usually difficult to fulfill (for a discussion of related numerical issues 
see for example |16ll37j ). Therefore, it would be useful to have a different approach. In principle, 
a possible way to determine the initial state |i/)o) in terms of the post-quench creation-annihilation 
operators Z, Z^ is to implement the following program: 

• for an arbitrary value of the parameter Q, find initially the relation between the ZF Z 
operators of the theory and the physical field operator (f>, i.e. — f{Q; Zq). 

• use the continuity of the field as boundary condition in the quench process Qq — > Q 

fiQo;ZQ„)^fiQ;ZQ) (9) 
for deriving the relation between the old and the new ZF operators 

Zq„ = F{Qo, Q; Zq) = f-\Qo; f{Q; Zq)) (10) 

• write the initial state \Qq), which is known in the pre-quench ZF basis (and is typically the 
ground state defined by Zq^I^Iq) ~ 0), in the new basis using the above relation. 

If this program can be realized, the time evolution of the initial state in the new basis can be 
computed easily. Going in more detail, the first step of this program consists in expanding the 
physical field operator as a series in the ZF operators using all of its form factors, i.e. the matrix 
elements of the field (j>{x) in the asymptotic states. The second step involves the inversion of 
this series; this might require an ingenious ansatz for the function F. The third step requires to 
deal with the most general expansion of a state in the post-quench basis and to determine the 
coefficients of \^lo) term by term from the equation F{Qo, Q; ZQ)|f7o) = 0. 

While the first step is essentially a re-expression of the body of information obtained by the form 
factors program, the other steps are in general highly nontrivial. In order to partially circumvent 
these difficulties, in the following we will exploit some general properties of integrable field theories. 
As previously said, in free theories the relation between the new and the old creation/annihilation 
operators is Bogoliubov-type, fixed by the condition of leaving invariant the (trivial) ZF algebra 
of these theories. Analogously, for generic integrable theories, the transformation between the 
pre and the post-quench ZF operators must respect the algebra. This leads us to investigate 
under which condition this requirement is satisfied. Of course this is quite an abstract point of 
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view: knowing that a certain transformation respects the algebra does not necessarily clarify the 
physical nature of the quench protocol. Nevertheless, it is surely important to understand what 
are the possible algebra-preserving transformations and if their form is restrictive enough to make 
prediction about the initial state. In Section 21 we integrate this analysis with a perturbative 
study of a typical integrable model, the Sinh-Gordon model, and its non-relativistic counter-part, 
i.e. the Lieb-Liniger model. 

3.1 Conditions required for transformations of the ZF algebra operators 

We are looking for transformations of the creation-annihilation operators Z'^ that respect the 
ZF algebra. We also require that the transformations respect the translational invariance of the 
theory, since we are considering only homogeneous systems, both before and after the quench. For 
this reason we will write the ZF algebra in momentum representation which, as we will see soon, 
serves better this requirement 

ZptZp^^ S[pi,p2)Zp^Zp^ (11) 
ZpiZl^_ = S{p2,Pi)Zl^Zp^ +S{pi -P2) , (12) 

along with the standard properties of the S'-matrix 

S{pi,P2y^ = S{pi,P2)* = S{p2,Pi) = S{-pi,-p2) . (13) 

Notice that in comparison with the form of the ZF algebra in the rapidity representation ([3]), we 
have redefinecjfl the operators as Zp = Z{p) = Z{9{p))/ E{p) since d{6{p)) ~ E{p)6{p). 

We focus our attention on infinitesimal transformations, assuming that a finite transformation 
can be build up by repetitive action of the infinitesimal ones. We also allow for infinitesimal 
changes of the S'-matrijQ. Since we demand that the transformations commute with the momentum 
operator, the new operator must carry the same momentum as the old one but not necessarily 
the same rapidity, as the quench may involve a change of the mass of particles (this is the reason 
why the momentum representation suits better our problem). Therefore both the transformed 
operator and the S'-matrix are in general expressed as 

Z; = + eWp (14) 
S'{pi,P2) = S{pi,p2) + eT{pi,p2) (15) 

where e is a small quantity, function of the infinitesimal change SQ of the quench parameter. In 
order to satisfy the ZF algebra, they must fulfill the conditions 

Wp.Zp, + Zp.Wp, ^ Tip,,p2)Zp,Zp, + S{pi,p2)iZp,Wp, + Wp,ZpJ (16) 
Wp.Zl + Zp.Wl = T{p2,pi)ZlZp, + S{p2,pi){ZlWp, + W^Zp,) (17) 

for all values of p,p' , along with the following conditions for T 

T{pi,P2)* = T{p2,pi) = T(-pi,-p2) = -T{pi,p2)S^'^{pi,p2) (18) 

coming from the unitarity of the S'-matrix. 

^This definition is also tailored to our purposes, since the energy of particles may change under a quantum 
quench and we would like to absorb all changes into the transformation of the operators. In addition it is consistent 
with the usual normalization of energy eigenstates in the free limit which will be useful later. 

''Note that this assumption may exclude the special case of free bosons. This is because for all integrable field 
theories except for free bosons, the 5-matrix at zero momentum is 5(0) = S{p,p) = —1. Therefore the transition 
from a free bosonic point of an integrable theory to another point that does not correspond to free bosons is 
always discontinuous as far as the S-matrix is concerned. In the following we sometimes make use of the property 
S(0) = —1 in which case we mention it explicitly. 
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The operator Wp can generally be written as an expansion in the operators 



OD / 71 rn \ n m 

-5 p + E * - E n n (19) 

The above conditions are then translated into a sequence of relations between the coefficients an,m 
of different orders. Below we construct and study several simple classes of solutions in which the 
above expansion terminates after a few terms. 



3.2 A first trial: linear Bogoliubov transformations 

Let us initially assume that W corresponds to a linear Bogoliubov transformation which, in the 
infinitesimal form, means Wp = apZ^_p. In this case it is easy to see that the previous conditions 
become 

S{-p,p') ^ S{p, -p')^Sip',p) 
ttp/a^p = S{p,p) = S{0) 
T{p,p')^0 

for all values oip,p'. The first of these equations implies that S{p^p'Y = li i-e. S{p,p') = ±1. We 
therefore arrive at the interesting result that the linear Bogoliubov transformation is a symmetry 
of the algebra only in the trivial case of free fields, bosons or fermions. Moreover, it is easy to 
show that any other linear combination of the operators is inconsistent with the general conditions 
(HH) and (dlD. 



3.3 Generators of S-matrix changes 

We will now construct a transformation that induces a non-zero change T in the S'-matrix and 
show that this transformation is unique, in the sense that any infinitesimal transformation that has 
the same effect must necessarily involve this one. First, observe that, since linear transformations 
are already excluded, the TZZ term in p6p can only be produced as a (5-function by-product of 
the commutation of higher order terms in W. More preciselly, W must be of 3rd order and must 
contain one operator and two Z's, so that the commutation of Wp with Zpi produces a ZZS 
term. Furthermore, from eq. p6p we see that the two Z operators in the residual ZZS term, which 
come originally from Wp^ must carry momenta p,p' (the same as Wp and Zpi) and therefore the 
Z'^ operator in Wp must carry momentum p' to ensure that Ws total momentum is p. Thus Wp 
is of the form Z^iZp/Zp. But this must be true for arbitrary p' , so W should necessarily consist of 
a linear combination of all such terms. All this leads to the ansatz 




Let us verify it explicitly by substituting into the required conditions (jl6p and p7p . the first of 
which gives 

T{pi,P2) = 5(pi,p2)(ap2.pi - api,P2) : (21) 

while the second 

T{p2,Pi) = S{p2,Pi){ap^p^ +api,p2) : (22) 

together with 

ap,q + a*p q = . (23) 

Remarkably all of these requirements are simultaneously satisfied as long as it holds the condition 
(|23p. i.e. if the coefficients ap^q are purely imaginary. Hence, from now on we set ap^q ~ iop q 
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where real functions. Notice that this solution also ensures that the S'-matrix remains 

unitary, as expressed by the conditions ((TH]) for T. 

Studying in more detail these transformations, one realizes that 



Z''^pZ p — ZpZp + ieZp 



'^{ap^g + a*pg)ZlZg 



L q 



Zp — ZpZp 



i.e. the conserved charges 



Qs = qs j dOe'^'zlZg 



(24) 



(25) 



remain invariant (unless the factors qs depend explicitly on the physical parameters whose in- 
finitesimal change leads to this transformation) . This allows us to easily derive the corresponding 
finite transformatiorlfl 

Z'p = ^{cxp [i j ip{s)ds^ }Zp where Jp(s) = ^ ap,,(s)ZtZ, (26) 
which changes the S'-matrix as 



S'{p,q) = cxp 



i / {aq.p{s) - ap,g(s)) ds 



S{p,q) 



(27) 



Even though we have constructed this infinitesimal transformation heuristically, it is easy to 
show that this transformation is the only one that changes the S'-matrix. Any other transformation 
that changes the S-matrix must necessarily be a linear combination of this one along with some 
other part that does not change it. Indeed if there was another transformation W' that also shifts 
S to the same S + eT then from ([T5)) and (|17p their difference W — W would not change the S- 
matrix. We can therefore decompose any transformation that respects the ZF algebra (infinitesimal 
or finite) into two parts, one of which is of the above form and performs the shift of the S-matrix 
to the desired value, while the rest leaves it invariant. In this way we have reduced the problem 
of finding the symmetries of the ZF algebra to the task of identifying those transformations which 
do not alter the S-matrix and which satisfy and ([T71) with T — 0, i.e. 



Wp.Zp, + Zp.Wp, = S{pi,p2){Zp^Wp, + WpM 
Wp.Zl + Zp.Wl = S(p2,pi)(Zt^W^,, Wl_Zp,) 



(28) 
(29) 



3.4 Other classes of transformations 

Having reduced the problem to identifying the transformations of the ZF operator which do not 
alter the S-matrix, we will now consider more general classes of symmetry transformations of the 
ZF algebra. Let us assume initially that W is simply a single product of Z, operators. In order 
to check the condition (j28p . we have to consider how Wp commutes with Zpr. when we swap Zpi 
with each of the operators in W one-by-one, this operation gives as output, for each of these terms, 
multiplicative S-matrix factors as well as additive (5- functions for each Z^, which are lower order 
products. One obvious way to satisfy (pS)) is then to choose Wp in such a way that (a) the overall 
S-matrix factor is simply equal to S{p,p') and (b) that the residual lower order terms vanish. 

Let us firstly focus our discussion on the point (a). We assume that Wp consists of n Z''- 
operators and m Z-operators (in some ordering that is not relevant for the moment), i.e. 

n ?n m n 

= n n ^i*^^ E - E = p • (30) 

1 = 1 j = l j = l 4=1 

^'Vl-..} denotes a path-ordering integration and s is a continuous parameter along some path. 
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Then we have 

(n 771 \ 

n ^(p'' n ^('^^ '^^') ^p'^p ' (31) 

1=1 3=1 ) 

where we use the symbol ~ to denote equahty for the highest order terms only (i.e. we ignore for 
now all residual lower order terms). To satisfy (|28p we then require that it holds the equatiorf^ 

n m 

Y[S{p\q,)l[S{r,,p') = S{p,p') . (32) 
i=i j=i 

In order to satisfy this relation for all p, p' and also independently of the specific functional form 
of the S'-matrix, we have to exploit its general properties. In particular, taking into account that 
S{p,p') = S^^{p',p) we see that if 

m = n + 1 and qi ~ rj for all i ~ j = l...n and r„+i = p (33) 

then eq. ([5^ becomes an identity. 

As for the second condition (j29p . we have 

(71 m \ 

\{S{q.,,p')\{S{p',rMzl,Wp, (34) 
1=1 j=i J 



and 



^pK ^ I n ^(^^ ' n ^(P' I K ^p ' (35) 

I j=l i=l 



and so we would similarly require 



l[S{q,,p')l[Sip',r,) ^l[Sir'^,p)l[Sip,ql) = S{p',p) . (36) 

i=l j=l j=l i=l 



Remarkably this condition is essentially identical to the one of eq. (|32|) . i.e. our solution ([33]) of 
([5^ automatically satisfies this one too. Hence, there exists a solution to these equations for 
arbitrarily high order n + m. 

However this is not the end of the story, since one has also to check the point (b), namely that 
the residual terms vanish. In order to ensure this condition for all p,p', instead of considering a 
single product pop . one has to look at a linear combination of such terms for all momenta qi and 
choose their coefficients so that the residual terms cancel each other. If this could not be realized, 
one would still have the option to introduce into Wp suitable lower order terms and cancel the 
residual terms order by order. In this way the coefficients of terms of order n depend on those of 
order n + 2 and we see that the construction of the transformation can be carried out recursively. 

Provided that all of these requirements are met, the resulting transformation is of the form 

/ „ \ 

Wp = a{p, {qi}) j ^li^qi ] + (suitable lower order terms) (37) 
{9.} Vi=i / 

Here we only report the lowest order members of this family of transformations. The first one 
corresponds to ?i = 1 and is the one we have found already in the previous section (in this case 
the residual terms result in a nonzero T, as we saw) 

Wp = i^^ap^qZ^ZyZp . (38) 

9 



®Of course this is not the only way to meet the condition II28II under II30I but, as the more detailed discussion 
presented in appendix^ shows, the other options lead, at the end, to the same form. 
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For n = 2 we find that the coefficient must be simply an imaginary constant 

Wp^iY, ZlZgZlZrZp . (39) 

g.r 

Our study started by assuming that Wp is a single product of Z, Z^ operators, a monomial 
(even though later we had to generalize our assumption by considering linear combinations of 
similar terms and lower order ones). However this is obviously not the only possibility. Another 
possibility is investigated in appendix |Bl which starts from a binomial and leads to the discovery 
of another interesting type of transformations 

Wp = Y, b,iSp,,Sp^,, - l)ZpZlzlg + ^ 6*(5p,,5p,_, - l)Z^,Z,Zp + 2h^pZlp = 
q q 

= bg{l - Sg,pS.g,p)Zlzl^Zp + KiSp,qSp.-g - \)Z-^Z^Zp - 2bpZip , (40) 
q q 

where bq has been chosen to satisfy b-q = bqSq^-q and we have assumed that the S'-matrix satisfies 
5(0) = Sp^p = —1. As already mentioned this excludes only the case of free bosons since for all 
other integrable models it is always true. Note that the last term cannot be absorbed by reordering 
the operators of the first one. 

Let us remark that one may continue in a similar way and construct other more complex classes 
of transformations. In particular, one may even consider the infinite series of products (jl9p which, 
unlike all cases presented above, do not give rise to expressions that terminate at finite order. The 
study of such transformations will be discussed elsewhere. 



3.5 Properties of the two simple classes of transformations 

In the previous sections we have found mainly two distinct classes of symmetries of the ZF algebra 
which led to eqs. (PU)) and PO)) . The first type of transformations (PU)) is the one that generates 
a change in the S'-matrix. However this first type does not change the ground state of the theory 
since Z'p = Zp + eWp annihilates the same vacuum as Zp. The reason is that this transformation, 
as well as the second member of the same class ppp . contain always one more Z-operator than 
Z^'s. Finally, let's notice that it does not reduce to the Bogoliubov transformation in the free 
limit S = ±1 since it does not depend explicitly on S. 

The second type of transformation PO)) has three important properties. Firstly, it does not 
change the 5-matrix. Secondly, in the free limit where S ±1, its nonlinear terms (first and 
second) vanish, leaving only the linear term (last) Z^_p, which corresponds to the Bogoliubov 
transformation. Thirdly and most importantly, it changes the ground state, since the first and 
last terms in (|40p contain one more operator than Z's which means that the new annihilation 
operator does not annihilate the old ground state. In particular, as we will show next, the in- 
finitesimal change in the ground state can be described as creation of a pair of excitations with 
opposite momenta. 

Indeed, if we denote by \rt) the ground state corresponding to the pre-quench operator Z, by 
definition this state satisfies 

Zk\n) = (41) 

for all k. This condition, expressed in the basis of the post-quench operator Z' (with corresponding 
ground state |fi')), reads 

(Z[. - eWk){l + eX)\n') = e{Z'^X ~ Wk)\n') = (42) 

where X is a suitable operator to be determined. For W given by P0|) we easily find by normal 
ordering that 

X ^ ^Y^pZ'lz'lp . (43) 
p 
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Figure 1: Diagrammatic representation of the three types of transformations in the context of quan- 
tum quenches. The linear Bogoliubov transformation that works only in free systems corresponds to the 
transformation of an old particle into a new antiparticle of opposite momentum. The two new types of 
transformations that work m the interacting case, convert the old particle into a new particle of the same 
momentum accompanied (or "dressed") by a particle- antiparticle (type I) or particle-particle (type II) pair 
with opposite momenta. 



Notice that this is the infinitesimal version of a squeezed coherent state. In fact we can go 
much further and show that for real bp any finite transformation generated by eq. (|40p "transforms" 
the initial ground state into a squeezed coherent stat^ll- To prove this it is sufficient to show that 
any state j'l') of the squeezed form 

I*) = M{K,) cxp K.Zlzi^^ (44) 

in the pre-quench basis preserves its squeezed form under the infinitesimal transformation (|40p . 

'^By "transformed" ground state, we mean the expansion of the ground state of the pre-quench operator Z into 
the basis of the post-quench operator Z' . 
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i.e. it is "transformed" into a squeezed state in the post-quench basis 



(45) 



If this is true then, since the finite transformation can be built up by successive application of 
infinitesimal ones and the initial ground state is "transformed" into a squeezed state (|43p which 
remains of this form after every infinitesimal step of this procedure, we conclude by induction that 
the transformation (PU|) indeed "generates" squeezed states. 

To prove this statement we can follow a path parallel to the corresponding free field calculation. 
For free bosons or fermions one should show that under an infinitesimal Bogoliubov transformation 



Z'p — Zp 



eapZ_p, 



{a-p = Gp) 



(46) 



any squeezed state of the form (|44|) remains squeezed as well. The easiest way to see this is to 
employ the following equivalent form of ([^^ 



cxp 



Z-qZq 



\n) 



(47) 



where is a known function of Kg 



In the above we have assumed that a^, Kq and are all real 
functions and all of the following are restricted to this case. The operator (Z^zl^ — Z^qZq) in 
the exponent remains invariant under the Bogoliubov transformation (j46p and therefore the only 
change comes from the ground state 



\n) = (^i-eY,aqZ'lz'lg^ \ny 



which can be absorbed in a shift of the coefficient A^ in the exponent 



I*) = exp 



^{Aq - eaq) [Z'lz'l^ - Z'_^Z'^ 



(48) 



(49) 



This is exactly what we wished to show. To verify that eq. (|47l) can also be written in the form 
dU]) one can normal order the squeezing operator exp[^^ kq(Z]^z'^_g — Z^qZq)\. Alternatively one 
may first decompose the exponential of the sum over momenta in (j47p into an infinite product of 
exponential^. Then observe that for each pair of opposite momentum modes, the commutation 
relations of the operators Zjzl^, Z-qZq and ^{Z^Zq + Z^_gZ-q±l) (+ for bosons / — for fermions) 
form a closed algebra {SU{1, 1) algebra for bosons / SU{2) algebra for fermions) and therefore we 
can rewrite the exponential as 

g2A,(Ztzl,-Z_,Z,) ^ g2A:,(A,)Ztzi^ g2L,(A,)Z,Z_, gA/,(A,)(ztz,+zi^Z_,±l) ^qq^ 



for some appropriate coefficients Kq, Lq, Mq that it is not necessary to determine here. Applying 
this operator to jfi), we directly derive the result mentioned above. 

Extending this computation to the general integrable case is almost straightforward. The first 
step is to check that even in this general case, squeezed states can be equivalently written in 
both forms and (|T7|) . This is true, since the commutation relations of the operators involved 
in the computation are not crucially different from the free fermionic ones. Indeed, if we define 
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Z\Zlq, Nq 



Z\Zq and Nq 



ZqZ\ 



-A^„ + 1 then 



^q^^■p\ ~ ^q.-pS-p,p{N^p 
[Y^lNp]^-{Sq,^p + Sq,p)Y^^ 



Np)+Sq,p{Np~ N^p) 



(51) 
(52) 
(53) 



°To avoid ordering problems in this product, first restrict the summation variable q to positive values only, as 
we are allowed to do. 
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and the only difference with free fermions is the factor S'_p.p in ([5^ . The next step is to check if 
the operator (Z^zl^ — Z-qZq) in the exponent of (j47p remains invariant under the transformation 
which turns out, after some algebra, to be true if hp is real. This completes the proof of our 
statement. 

Let us now discuss another property of this second type of transformations: the change induced 
to the conserved charges of the theory. As we have already seen in , the first type we studied 
leaves all conserved charges invariant. In the present case the transformation of Z^Zp turns out 
to be 

Z'lz'p - ZlZp - 2e(6pZtzlp + b;Z^.pZp) , (54) 

that is the conserved charges Qg do change but they remain quadratic in the Z operators. In 
particular the hamiltonian of the system which in momentum representation and according to our 
normalization of the ZF operators is iJ = EpZj^Zp is transformed as 



H' = H + eJ2 -gfK^P - "^EpihpZlzl^ + h^Z^pZp) 



(55) 



where dEp/de reflects the change in the dispersion relation induced by the transformation. Notice 
that this expression is reminiscent of the analogous one for Bogoliubov transformations in free 
theories. 

From a physical point of view and in particular from the perspective of quantum quenches, 
the action of transformations of these two types consists in aAYdressingaAZ the initial particle 
of the theory with a pair of newly created particles with opposite momenta, as depicted in figure 
[TJ Especially for the second type, if such a transformation describes a quantum quench, then the 
initial ground state is expressed as a squeezed coherent state in the post-quench basis, at least for 
real hp. Finally, the perturbation introduced in the hamiltonian after a small quantum quench of 
this type would be of quadratic form in the pre-quench ZF operators. 



4 Examples from physical theories 

In this section we will present some first examples of transformations of ZF operators in the context 
of two important integrable models, the Sinh-Gordon and the Lieb-Liniger. The first is one of 
the simplest and best-studied relativistic integrable models consisting of a single type of particles, 
while the second describes a system of non- relativistic interacting bosons in (l+l)d and models 
experimental cold atom set-ups. As it has been recently shown in |38| . the two models are closely 
related to each other, since the Lieb-Liniger can be obtained as a suitable non-relativistic limit of 
the Sinh-Gordon model. 

In the Sinh-Gordon model, an example of infinitesimal transformation would correspond to 
a small quench of the mass m or the coupling constant g. Starting from arbitrary initial val- 
ues, such a quench would introduce into the hamiltonian perturbations (J dx cosh g(j)[x) and 
J dx (2 co&h- g(j){x) + g(l)(x) smhgcjjix)) respectively) that correspond to an infinite series in terms 
of the Z, Z'^ operators, as can be seen from their form factors 131]. An exception to this rule is 
when the initial model lies on a free point, 5 = 0. In this special case where a small g is abruptly 
switched on, each term of the expansion is smaller than the previous one by an amount of the 
order of g^, due to the same property held for the form factors ^"2^1+1 ^^'^ physical field. Thus 
the derivation of the corresponding infinitesimal transformation is simpler and can be done by 
means of perturbation theory, which is what we do in the next section. In the Lieb-Liniger model 
on the other hand, we use the known relation [lU] between the physical field operator and the ZF 
operators to find the infinitesimal transformation of the latter, again for the case when the inter- 
action changes from zero to a small value. Using the non-relativistic limit mentioned above [3H1 
we can verify the consistency of the results for the Sinh-Gordon and Lieb-Liniger models. 

Note that in both cases the transformation refers to a small change from a free bosonic to an 
interacting point of the theory and, according to a comment we made in a footnote of 13.11 the 
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S'-matrix at zero momentiiin S{0) is non-analytic (in fact discontinuous) in the coupling constant 
at such points. Therefore we anticipate this non-analyticity to become evident in our results and 
indeed it does as we will see below. 



4.1 The Sinh-Gordon model 

The Sinh-Gordon model is a relativistic field theory in (l+l)d defined by the hamiltonian 

where (f> = (/)(a;, t) is a real scalar field, to is a mass scale and c is the speed of light. In this 
integrable field theory there is only one type of particle with physical (renormalized) mass M 
given by 

oSinaTT . . 

M=m^ (57) 

an 

where a is the dimensionless renormalized coupling constant 

a = (58) 

Particle scattering is fully determined by the two-particle S'-matrix given by 

sinhg - i sin an 

Ssi,.G{0,a) = (59 

smh A- 1 sm an 

where is the rapidity difference between the particles. 

To calculate the ZF operators from first order perturbation theory we consider the (j)^ model 

with coupling constant A ~ rn^c^g^ (we set c = 1 from now on). We first define the auxilliary 
operators 

B\{k) ^ Vt+A\k)Vt\ , Bl(fc) = 17_yl^(fc)0L (61) 
where f2± are the following evolution operators (M0ller operators) 

17± = lim e-~-S-TdtH.r.At) /g2) 

T~i±oa ^ ' 

As known from the general scattering theory, the operators i?^(fc), _Bl(fc) when acting on the 
vacuum state of the interacting theory jfi), create "in" and "out" states respectively. If we consider 
the interaction hamiltonian Hint as normal-ordered, we have 

Bi(fc) = At(fc)-^A ( dt ( dx[:(j}\x,t) :,At(fc)] = 

= A\k)-i- dt dx ■.^^ix,t): (63) 



and expanding in terms of the free boson creation/annihilation operators A{k), A^ (k) 
Blik) = At(fc)- 

g-iB.t r dkidk2dk32nS{k + '£^ h) 



dt- 



A(fci)A(fc2)A(fc3)e-'(-^'=i+^'=2+^'=3)*+ 



3!;^^ J {2nY ^2^ Ek,Ek,Ek, 

+ 3At(-fci)A(fc2)A(fc3)e-''-^''i+^''2+^'=3)* + 3ylt(-fci)At(-fc2)A(fc3)e~*(-^'"i--^'=2+^^3)* + 



4- At(-fci)At(-fc2)At(-fc3)e-'(-^'=i-^''2-^'=3)* 



(64) 
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In scattering theory the ^-integration is understood under the "adiabatic switching" prescription 
which means to introduce an e~'^'*' factor into the integrand. According to the identity 



■ 

dt e^^l*!-*"* = 

±00 wT^e 



(65) 



we then find 
Bl{k)=A\k) 



X 1 



dkidk2dk^2'!i5{k + ki) 



A(fci)A(fc2)A(fc3) 

Ek + + + Ek3 



3!^/2^ 

, At(-fci)^(fc2)A(fc3) 
} -r o :— -1 

Ek — Eki + Ekr, + -Efc., Ek — Eki — Ek2 + Ek^ ±ie Ek — Ek 



{2Tif^2^Ek,Ek,Ek, 

At(~fci)v4t(-fc2)A(fc3) , A\-kMH-k2)A^{-h) 



E, 



k2 



Ek 



(66) 



Notice that we kept the ±ie shift only in the 3rd term since this is the only one that has a 
singularity (at fci = — fc, ^2 — —k^ or k2 = —k,ki = —k^). Using the formal identity 



lini 



1 



i.o+ u! zt ie 



= V 



=F ITT (5(cj) 



we can rewrite the above as 



Bi{k)=B\k)T 



iX f dq 



1 



27r \kEq - qE,, 



rA\k)A\q)Aiq) 



(67) 



(68) 



where 
Bt(fc) = A^k) 



X 1 



, At(-fci)A(fc2)A(fc3) 

Ek — Eki + Ek2 + Ek3 



dkidk2dk^2'!i5{k + ki) 
{2tiY .JWWk.Ek^Ek, 



A{k^)A{k2)A{k^) 
Ek + Eki + Ek2 + Ek3 



3-p.V, 



At(-fci)At(-fc2)^(fc3) , At(-fci)At(-fc2)At(-fc3) 



Ek — Eki — Ek2 + Ek.^ 



Ek — Eki — Ek2 — Ek^ 



(69) 



From their definition ((6T|) . each of the two operators i?j-(fc) satisfy the same standard commu- 
tation relations as the free operators A{k),A''{k). Now let us define the operator 



Zt(fc) = B^k) 



iX f dq 



1 



2tt \ kEq - qEk 



iX f dq 
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1 



A\k)A\q)A{q) 
1 



27r \kEg-qEk \kEq-qEk 



A\k)A\q)A{q) 



(70) 



Since this is of the form (j20p . we automatically know that Z\k) and Z(k) satisfy the ZF algebra 
with non-trivial S'-matrix given by 



S{k,q) = l- 



iX 



A{kE, - qEk) 

which is indeed the correct first order perturbation to the 5-matrix ((59l 



(71) 



5sli-G \ 



9, a) ^l-cg"' 



4 sinh 9 



(72) 



Notice the infrared singularity in the coefficients of ([70)) and ([7^ when the momentum difference 
k — q tends to zero. This reflects the non-analyticity of S'(O) as g — > for which we talked in the 
introduction of this section. 
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Next we consider the states created by the action of Z^{k) on the perturbed vacuum \rt) which 
by a calculation similar to the ones above, turns out to be 

For the one-particle states it can be immediately seen that Z^(fc)|il) = i?j_(fc)|0) always up to 
first order in A. However, in order to verify that Z'^ {k) plays the right role in creating in and out 
scattering states, we should check the two-particle states Zt(fci)Zt(fc2)|51). By normal ordering 
we find 

z\k,)z^ik.m = Biik.)Biik.m - ^ (^^^^^^ T i.^^Jfe^^^J A^ikM^km 

(74) 

Observing that k/ Ek is a monotonically increasing function of fc, we can easily see that, if fci > ^2 
then 

Zt(fci)Zt(fc2)|r!) = Bl{ki)Bl{k2)\n) (75) 
i.e. it defines an in state, while if ki < k2 then 

Z^iki)Z\k2)\n) = Bliki)Bl{k2)\n) (76) 

i.e. it defines an out state. 

This example, apart from demonstrating how ZF operators emerge from standard perturbation 
theory of scattering, it also illustrates the concepts developed before and in particular the role of 
the first class of transformations that we derived abstractly. We close our presentation of 
physical examples with the Lieb-Liniger model. We will also verify the consistency of our results 
for the two models, under the double non-relativistic limit that reduces the former to the latter. 



4.2 The Lieb-Liniger model 

The Lieb-Liniger model describes a (l+l)d system of non-relativistic bosons interacting with each 
other with a ^-function potential. Its hamiltonian in second quantized form is 

H = J dx (^^d^^\x)d^'i'{x) + X^''{x)^''{x)^{x)^{x)J , (77) 

where A now is the interaction strength. The ground state energy for a system of N bosons as 
well as its thermodynamics can be exactly worked out by means of Bethe ansatz [41] . The exact 
solution expresses the energy of the ground state and the excitation spectrum in terms of the 
dimensionless coupling constant 7 = A/p where p = N/L is the density of bosons with L ^ 00. 

In this model the relation between the bosonic field operator ^'(a:;) and the ZF operators R\{k) 
that diagonalize the hamiltonian for L — > 00 has been already found using the inverse scattering 
method gU] 

CO N , W J7, 

*(^) = E / n ^ n ^ 9n{{p}, {k}; x) Rlip,) ■ ■ ■ Rl{p^)R,{k^) ■ ■ ■ R,ik,)R,{ko) (78) 

AT n " ; 1 n 



7V=0 i=l j=0 

where 



9N{{p},{k};x) = —j^j— ; ; — . (79) 

11^=1 (Pj - fcj - - kj-i - «e) 

Indeed it can be shown that the R, i?^ operators diagonalize the hamiltonian and satisfy the ZF 
algebra 

[H,Ri{q)] = q'Riiq), (80) 

Rx{q)Rx{q') = Sx(q' -q)Rx{q')Rx{q), (81) 

Rx{q)Riiq') = Sx{q-q')Rliq')Rx{q)+27rSiq~q') (82) 
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where the iS-matrix is 

Sx{q) = . (83) 

Let us consider the infinitesimal transformation from the free bosonic point A = to a smaU 
value A. From ((75|) we hav^ 



R,{k) = Roik) + X T V + 0{X') (84) 



dqdq' Rl{q + q' -k)Ro{q)Ro{q') , ^^^^^^ 
(27r)2 {q - k - ie){q' - k - ie) 

The S'-matrix is no longer unit but becomes instead 

S^{p) = l-X^ + 0{X') (85) 

and so, according to our previous findings, we expect that the transformation contains the gener- 
ator of S'-matrix shifts ([^0]) with coefficients ak,q = i/{q— k). Indeed, using the identity (|67|) we 
recognize that part of the infinitesimal transformation (|M)) has exactly the form of ([20)1 with the 
right coefficient 

'^'^ ' -Rl{q)Roiq)Roik) , (86) 



277 q — k 

while the remaining part does not affect the S'-matrix. Once again, notice the infrared singularity 
in the coefficient of the above expression for q ~ k. 

Lastly, we mention that the infinitesimal transformations (|70p and (j84p derived for the Sinh- 
Gordon and Lieb-Liniger model respectively, are consistent with each other under the double 
non-relativistic limit c — >■ oo, g 0, gc : const., that leads from the former to the latter model. 
Following [3S| we substitute the field in as 

</.(x, t) = f V(x, i)e-™^'* + ^\x, i)e+™='*) (87) 

and keep only the non-oscillating terms, rewriting first all expressions with their c dependence 
explicit and taking into account that in the non-relativistic limit Ek = mc? + k"^ /2m + .... After 
some algebra we verify that ([70)1 reduces to ([M]) . 



5 Conclusions 

In this article we have investigated how the initial state of a quantum quench in an integrable 
model can be expressed, from first principles, in terms of ZF operators, without relying on the usual 
mapping to slab geometry and the associated boundary renormalization group arguments [1111121 
[23]. We show that this result can be achieved by deriving the relation between the pre-quench 
and the post-quench operators on the condition that such a relation respects the Zamolodchikov- 
Faddeev (ZF) algebra satisfied by integrable models. 

Under the conditions that such transformations must satisfy the ZF algebra at the infinitesimal 
order, we have initially showed that the usual linear Bogoliubov transformations do not respect 
the ZF algebra, apart from the trivial cases of free bosons or fermions, a result that holds generally 
for finite transformations too. 

We have then identified two important classes of transformations. The first class changes 
the S-matrix of the theory but preserves its ground state as well as its conserved charges. We 
also argued that any infinitesimal transformation can be decomposed into a part that induces 
the 5-niatrix shift and the rest that does not alter the S-matrix. The second class belongs to 

^Note that, unlike in a relativistic free field theory where the creation-annihilation operators are linear combi- 
nations of the field and its conjugate momentum tt, in a non-relativistic free field theory the creation-annihilation 
operators are the bosonic field itself Ro{k) = J dx e~*'''^^'(a;). Also the conjugate momentum is 11 = i'tt and does 
not appear in the hamiltonian. 
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the latter subset of transformations, which can be regarded as a generahzation of the Bogohubov 
transformations for interacting theories, since it reduces to the usual Bogoliubov transformations 
whenever the integrable model reaches a free bosonic or fermionic point. Like in the free case, the 
ground state of the system becomes a squeezed state when expressed in the transformed ZF basis 
under such a generalized Bogoliubov transformation, at least when its coefhcients are real. We 
have also showed that the change in the hamiltonian (and in the other conserved charges of the 
theory) is of the same form as that of the Bogoliubov transformations. The net effect of this type 
of transformations is to "dress" the initial quasiparticle with pairs of new particles with momenta 
opposite to each other. 

We have also outlined how one could proceed further in this program to identify the transfor- 
mations which preserve the ZF algebra, in particular pointing out the existence of transformations 
of higher complexity characterized by the fact that, even for infinitesimal quenches, they are asso- 
ciated to an infinite series of terms given by products of the initial creation/annihilation operators. 
In the quantum quench perspective, this means that even a small quench of the physical param- 
eters of an integrable model may result in an infinite series which links the pre-quench and the 
post-quench operators. In this case, the calculation of the initial state made on first principles is 
rather difficult, unless a truncation or resummation of the series can be established on the grounds 
of a different argument. In such a case, for instance, it may be possible to reorganize the terms of 
the series based on a small-density expansion, following the concepts developed in [5S]. We hope 
that our work on ZF algebra transformations will stimulate further investigation of their structure 
and properties, both from a pure and an applied point of view. 

Lastly we exemplified our approach in the context of the Sinh-Gordon and Lieb-Liniger model. 
We restricted ourselves to perturbations about the free bosonic point of these models since in 
this case the transformations can be found relatively easily and contain only up to cubic terms 
in the ZF operators. We expect that analogous simplification occurs near free points of other 
integrable models too and it would be interesting to explore some physical realization of such 
quench processes. Regarding the Sinh-Gordon model, an application of our results to the corre- 
sponding quantum quench problem of an abrupt switch-on of the interaction, would give results 
comparable with earlier work |42| . For the Lieb-Liniger model instead, further manipulation is 
required, mainly due to the fact that the ground state is not the empty vacuum but contains 
a large number of particles proportional to the size of the system. This issue is discussed and 
a numerical approximation is developed in |37) . A recent numerical study of a special quantum 
quench in the Lieb-Liniger model is [17| . 
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A Squeezed states in free quantum field theories 

In this Appendix we show that the squeezed coherent states in a quantum quench of free theories 
are a consequence of the Bogohubov transformation of their operators. 

Bosonic theory. Let us consider firstly the quench in a bosonic theory with hamihonian [T^] 



H 



[n^ + {V(p) + mQip^l dx 



This system can be diagonahzed in momentum space 

H = 



qO . aO-\ AO . 



1 



with the ground state j^I'o) identified by the condition 

^°l*o) = . 

Imagine now that, after having prepared the system in its ground state, we quench the mass 
mo — >■ m. The relation between the pre-quench ladder operators {A^,A^^) and the post-quench 
ones (Afc, A^) is a Bogoliubov trasformation 



Ak 



CkAl + dkA% 



Al = CkAk - dkAl^, 



Ife - Cfe^fe 

ifc ~ ' 
where the coefficients are given by 

^k + 



A\ = + rf^^o ^ 

A\^ ^CkA\-dkA.k. 



Ck 



2^Jflk^ 



dk = 



^k 



Substituting the expression of A^ from the Bogoliubov transformation into eq. 
in terms of the new operators, the initial state satisfies the condition 



CkAk - dkA''_, 



l^n) = 



we see that. 



(89) 



whose solution is given in terms of a squeezed coherent state 

|*o) = A/'exp 

where 



A', 



hos 



oson 

ik)AlAl^dk 

nl - Qk 



|0), 



nl + nk 



(90) 



(91) 



This quantity can be written in a suggestive way by introducing the rapidities of the particle 
relative to the initial and final situations, i.e. 



il*^ = mo cosh 6*0 
ft = m cosh 9 



k = mo sinh 9q 
k = m sinh 9 



(92) 



19 



From the equality of the initial and final momenta, we have the relation which links the two 
rapidities 

Too smh6 

TOosmhwo = msmhy ^ = . (93) 

TO smh Oq 

and therefore, the amplitude Kboson{k) of eq. (|9ip can be neatly written as 

Too cosh 00 - "1 cosh 6^ ^ cosh 6*0 - cosh 6* 
rriQ cosh Oq + m cosh ^ cosh do + cosh 
sinh 6 cosh Oq — sinh Oq cosh 6* sinh(0 — 6*0 ) 



sinh 9 cosh 9q + sinh Oq cosh ^^ sinh(6' + 6'o) 

Fermionic theory. One can easily work out the Bogoliubov transformation relative to the 
quench of the mass of a free fermionic system [23] ■ Consider, in particular, a free Majorana fermion 
in (1+1) dimensions, with the mode expansion of the two components of this field given by 

/ + 00 
dp [a(p)A(p)e~*^*+*f" +a(p)At(p)e*^*-*f"] 
-OO 

/ + 00 
dp [p{p)A{p)e-'^'+'P- + /3(p)At(p)e'^*-'P-] 
-OO 

where 



u ^E+p uj ^/E + p 



27rV2 E ^ 27rV2 E 

with w = exp(z7r/4). At t = 0, i.e. at the instant of the quench, we can extract the Fourier mode 
of each component of the Majorana field ipii.^^ 0) = / dpipi{p)e'^'"^ , given by 

Mp) ^ a(p)A(p) +S(-p)At(-p) 

Mv) - mm^^{-v)A\~p) 

Suppose now that the mass of the field is changed from toq to to at i = and let's denote by 
(Ao(p): ^o(p)) ^^'^ the sets of oscillators before and after the quench. The proper 

boundary condition associated to such a situation is the continuity of the field components before 
and after the quench, i.e. ip'^{x,t = 0) = ipi{x,t = 0), which implies i/Jiip) ~ i^iip)- This gives 
rise to the Bogoliubov transformation between the two sets of oscillators 

Aoip) = u{p)A{p) + iv{p)A\-p) 
Alip) = u{p)A^p) - iv{p)A{^p) 



where 



1 



uip) = ^ [ViEo+p){E+p) + ^iEo-p)iE-p) 
1 



vip) = ^ [^iEo^p){E + p)- ^{Eo+p){E-p) 

Notice that these functions satisfy the relations u{p) = u{—p) and v{p) = —v{—p) together with 
u^{p) + v^ip) = EojE^ which refers to the normalization of the respective set of oscillators. 

With the same procedure used in the bosonic case, it is easy to see that the boundary state 
corresponding to this quench can be written as 

\B) = exp (J dp Kf,r^,on{v)A\p)A\-p)\^ |0) , 
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where 

As in the bosonic case, this quantity can be expressed in a more concise form by introducing the 
rapidities of the particle before and after the quench, i.e. 

Substituting these expressions in (pSj) . we get 

sinh I 



In conclusion, the squeezed coherent form of the initial state in free theories comes from the fact 
that, in any quantum quench of these systems, the creation-annihilation operators before and after 
it, are related by a Bogoliubov transformation. And this in turn is a consequence of the canonical 
commutation/anticommutation relations satisfied by these fieldj^. 

Finally, notice that both the bosonic and fermionic amplitudes K{p) do not satisfy, in general, 
a unitarity equation, in contrast to the amplitudes of squeezed coherent states in the purely 
boundary integrable theories studied by Ghoshal and Zamolodchikov |36| . The reason why this 
condition is requested for purely boundary field theory but not for an arbitrary quench process 
is quite easy to understand. In purely boundary theory, all the degrees of freedom for t < are 
completely frozen: the hard-wall boundary does not allow any transmission process through the 
boundary and this ends up in the unitarity condition. However, for quantum quenches in free 
theories, there are degrees of freedom also for t < 0, which are related to the ones at t > just by 
the Bogoliubov transformations. Said differently, the boundary condition in free theories allows 
for transmission, as it is shown in figure The degrees of freedom present on both sides of 
the boundary prevent the bosonic and fermionic amplitudes (j9ip and (|95p to satisfy a unitarity 
equation. Notice that the only case where they satisfy a unitarity equation is when we freeze the 
degrees of freedom before the quench by taking the limit rriQ — >■ oo: in this limit, in fact, we have a 
transmission-less boundary which implements in both theories the Dirichelet boundary condition, 
1 and Kfermion = * tanh I respectively. 




Figure 2: The boundary conditon for a mass quench allows for transmission. Notice that, for free 
theories, there is no particle production at the boundary and the transmitted particle has always the 
same momentum as the incoming one, since momentum is conserved at the boundary. However, 
since the mass is different at x < and x > 0, the rapidity changes from 6 to 9' . 



^"To be precise there can be exceptions to this rule since it is possible to construct generalized Bogoliubov 
transformations which satisfy the CCR/CAR but are nonlinear 1431145 1 . These however correspond to non-quadratic 
hamiltonians which, even though they can be reduced to free ones, they are uncommon in physically interesting 
cases. 
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B Derivation of two classes of generators of ZF algebra trans- 
formations 



In this appendix we present the derivation of the two simple classes of infinitesimal transformations 
of the ZF algebra that we introduced in the main text. As already mentioned we can always write 
W a.s a linear combination of products of Z and Z^ . The first class arises when we consider W to 
be a single such product 

Wp=l[z,l[z} with p = Y^p.^-J2p^ (97) 

For brevity we use the notation Zi = Zp. when there is no confusion about the meaning of the 
indices, and also the symbol '~' that we defined in paragraph [2131 Then we find 

W1Z2 ~ \12Z2W1 (98) 

where A12 = Yl Si2 Y\ S2j (note the implicit dependence on pi through the momentum condition 
in ((97)) ). Similarly we have 

WiZl ~ \l2ZlW1 (99) 
From (fro| and pT)) we have the conditions 

WiZ2 + ZiW2':^Sx2{Z2Wi + W2Zi) (100) 
WiZl + ZiWl ~ S2i{zlWi + wlZi) (101) 

where we ignored T12 since the corresponding terms are not of highest order (unless is a single 
operator instead of a product, which is the linear case that has already been excluded as not 
fulfilling the conditions). In the main text ea. p2|) . we considered the simple choice A12 = ^12 
which automatically satisfies the above condition but may not be the only possibility. Indeed 
poop is in general a weaker condition than this choice. Strictly speaking the two conditions above 
require 

(A12 - Si2)Z2Wi + (A*i - 5i2)W^2^i ^ (102) 
(AJ2 - S2i)zlWi + (A21 - S2i)WlZi ~ (103) 

which means, either A12 = 5i2 as before, or that W2Z1 and Z2W1 contain exactly the same 
operators and similarly for W^Zi and Z^Wi. In this last case Wp must contain Zp i.e. Wp = ZpWo 
where Wq has momentum zero and, as it can be easily deduced from the above equations, it 
must be independent of p and contain the same operators as its hermitian conjugate Wq, that is 
Wo -^Ylzlz,. But this is again exactly the same case considered before ([55)) and corresponding 

to A12 = >S'l2. 

Another possibility arises when we consider Wp to be a linear combination of two products of 
operators of the previous form 

Wp = a+,pVF+,p + a_,pVF_,p (104) 

The conditions are then 

a+iW+iZ2 + a+2ZiW+2 + a^iW-iZ2 + a_2^iW^-2 — 

Si2ia+iZ2W+i + a+2W+2Zi + a_iZ2VK_i + a_2VF_2^i) (105) 

a+iW+iZl + a*^^ZiWl2 + a^iW-iZl + ai^ZiWl2 ^ 

S2i{a+iZlw+i + al2Wl2Zi + a^iZlW-i + a*_2Wl2Zi) (106) 
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or 



(A+12 - Si2)a+iZ2W+i + {X*^2i - Si2)a+2W+2Zi + 

(A-12 - 5*12) a- 1^21^-1 + {X*_2i - Si2)a-2W-2Zi ~ (107) 
(A;i2 - S2i)a+iZlw+i + (A+21 - ^21X2^^+2^1 + 

(AI12 - S2i)a^iZlW-i + (A_2i - 52i)al2W^-2^i - (108) 

Like before one obvious choice is A+12 = A_i2 = S'12 which is trivial since each of VF+ and W- 
must fall in the earlier discussed case, but unlike before there is an alternative that is not equally 
trivial. If Z2W+1 and W-2Z1 contain the same operators and the same holds for Z^W+i and 
W^_2Zi (and similarly for the other pairs) then W+^p and W-^p can be written in the form 

W+,p = ZpWo and W-,p = ZpW^ (109) 
where Wq has momentum zero and must be independent of p. If we define 

ZpWo^fipWoZp (110) 
then A+12 = A'2'5'i2, A_i2 = H2S12 and the above conditions are equivalent to the set of equations 

a-2 ^J■2-^ 



a-i fH - 1 

a+2 ^ A^2 - 1 
a+i fi'l - 1 

«-2 _ M2 - 1 
a+i A^i - 1 

a+2 _ M2 - 1 
a_i ^1-1 

which have to be valid for any choice of pi,p2- The solution is 



(111) 
(112) 
(113) 
(114) 



a- 



,p = -b*{iip - 1) (116) 



for some constant h. 

A simple choice for Wq is the product of two operators that create a pair of particles with 
opposite momenta Z^zl^. Obviously a sum over all such pairs can be used, leading to the 
following expression for Wp 

Wp^Y. hiSp^iSp.-q - ^)ZvZ\z\ + ^ &;(5p,g5p,_, - \)Z_^Z^Zp (117) 

9 9 

where by symmetry of the sums under g — > —qwc find that bq can be chosen to satisfy h^q ~ 
hqSq-q. Substituting into the general conditions (|28p and (23) we find out that an additional lower 
order term proportional to is necessary in Wp in order to satisfy the full-form conditions. 
Overall the tranformation is 

Wp^Y. KiSp,qSp,-q - l)ZpZlZl^ + Y K^Sp.qSp^-q - \)Z^qZqZp + 2h-pZ\ (118) 

1 q 

Note that the last term cannot be absorbed by reordering the operators of the first one. 

Other choices for Wq are still possible and may lead to more complex families of transforma- 
tions. One may also consider linear combinations of more than two products and continue in a 
similar fashion. 
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